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Abstract 

Different metfiods for treating tlie results of higher-order perturbative 
QCD calculations of the decay width of the Standard Model Higgs boson 
into bottom quarks are discussed. Special attention is paid to the anal- 
ysis of the Mh dependence of the decay width r(H — > bb) in the cases 
when the mass of b-quark is defined as the running parameter in the 
MS-scheme and as the quark pole mass. The relation between running 
and pole masses is taken into account in the order a^-approximation. 
Some special features of applications of Analytical Perturbation Theory 
(APT) are commented. 



1 Introduction 

The study of the Higgs boson decay width into bottom quarks is rather impor- 
tant for calculations of the branching ratios of this important ingredient of the 
Standard Model and its various extensions. Current LEP and Tevatron fits of 
the Standard Model parameters yield the value for the Higgs boson mass around 
Mh = 76^11 GeV C.L. 68%. With the direct LEP search limit Mr > 115 GeV 
the fits provide at C.L. 95% Mh < 182 GeV. In the case, if the scalar Higgs 
particle has the mass in the region 115 GeV < Mh < 2Mw, its decay width 
to bottom quarks r(H — > bb) dominates over other channels. In particular, it 
is determining the branching ratio of H ^ 77 process which is considered to be 
one of the most promising channel for searches of Higgs particles at LHC in the 
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mass region specified above (for reviews see [1]— [4]). There is also a possibil- 
ity that the signal for H ^ bb process may be seen at Tevatron [5] through 
WH- and ZH- channels and at CMS-TOTEM [6] or/and FP420 [7] experimen- 
tal proposals at LHC, aimed at searches of central exclusive H production, as 
discussed from theoretical point of view, e.g., in Rcfs. [8,9]. The mentioned 
experimentally-oriented motivation is pushing ahead the intention to study in 
more detail the dominant theoretical effects to r(H — > bb) in the region of rel- 
atively light Higgs boson. These effects are related to high-order perturbative 
QCD predictions with their intrinsic uncertainties. Moreover, the comparison of 
various representations for r(H ^ bb) = T^^i,]^ is rather important for planning 
the experimental program of high energy linear e+e~-colliders for measuring 
Higgs boson couplings [4]. 



2 QCD expressions for Fjjbb 

Let us first consider QCD theoretical prediction for T^^^ expressed in terms of 
running b-quark mass and the QCD coupling constant in the MS-scheme as 



bmb(MH) 



i Hbb - i o 
™b 



l-f^ATiaUMn) 



i>l 



(1) 



Here Fq = Sv^/SttGfMh m^, m^, and Mh are the pole b-quark and Higgs boson 
masses, as(MH) = as{Mii)/iT and mb(MH) are the QCD running parameters, 
defined in the MS-scheme. The coefficients AFj are known analytically up to 
4-th order correction of perturbation theory [10]. They consist of the posi- 
tive contributions dp, calculated directly in the Euclidean region, and from the 
proportional to tt^ kinematic effects, which appear as the result of analytical 
continuation from the Euclidean space-like to the Minkowskian time-like region. 
This TT^-term arises first in Eq. (1) at the o^-correction [11]. Its coefficient was 
corrected later in [12], [13], but the kinematic vr^-term remained unaffected. 
Using the notations of Ref. [14] one can write down the following relations: 

AFi = df = ^; (2) 

AF2 = df -7o(/3o + 27o)vrV3; (3) 

AF3 = df - [df (/3o + 7o)(/3o + 270) + /3i7o + 271 (/3o + 27o)]7rV3 ; (4) 

AF4 = df - [df (/3o + 7o)(3/3o + 270) + df /3i(5/3o + 670 )/2 

+ 4df 7i(/3o + 70) + P270 + 271 + 71) + 72(3/3o + 47o)]vrV3 

+ 7o(/3o + 7o)(^o + 27o)(3/3o + 27o)7rV30, (5) 
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where the nf-dependence of df {i > 2) read 
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24 ~ 3 



51.567 - 1.907 nf « 42.032 (nf = 5) ; 



(6) 
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5015 
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432 



K 9470.8 - 1454.3 nf + 54.783 nf - 0.45374 nf «i 3512.2 (nf = 5) . (8) 

The term of Eq. (2) was evaluated in Ref. [11]. It is in agreement with the 
expressed in other ways results of previous studies, performed in [15-17]. The 
second coefficient was corrected in [12, 13]. The result of Eq. (7) was obtained 
in Ref. [18]. The exact value of the Euclidean coefficient of Eq. (8), analytically 
calculated in [10], turned out to be in reasonable agreement with the estimates, 
obtained within the used in Ref. [14] a variant of the effective charge approach 
(ECH) and the principle of minimal sensitivity (PMS) approach. The variant 
of the two approaches was developed in Ref. [19]. The coefficients /3j and 7^ 
enter the expansions of the QCD renormalization group (RG) /3-function and 
anomalous dimension of mass function 7^. The QCD /3-function can be defined 
as 



din ^j? 



= /3(as) 

= -/3o al - A a3 - /32 4 - ps ^ - /34 + 0(aj) 



(9) 
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Its MS-scheme expressions were calculated analytically up to 4-loop correc- 
tions [20], confirmed recently in the work [21]. We present here the results 
of analytical evaluation of the coefficients of Eq. (9) (Pq and Pi are scheme- 
independent) in the MS-scheme, supplemented by the numerical expressions, 
related to Uf = 5 number of active flavours: 
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/ 50065 6472 , 



n| + 



1093 



nf 



162 81 ^ J ' 729 
114.23 - 27.134 Uf -M.5824nf-F 0.0059nf «i 18.852 (uf = 5) . (13) 

Throughout this work we fix Uf = 5 and neglect the contribution of lighter four 
quarks to the relation between running mass rnb(MH) and pole (or on-shell) mass 
mb in Eq. (1) (more details will be given below). This is done for self-consistency 
of further analysis. Indeed, the similar contributions to the coefficient function of 
Eq. (1) are still unknown and are expected to be small. Notice an interesting fact: 
the growth of the coefficients of /^-function at Uf = 5 is starting to manifest itself 
from the four-loop only (on the contrary to the case with Uf = 3 when the values 
of the coefficients of perturbative series for the /^-function are monotonically 
increasing from the one- loop order). As to the five- loop coefficient P^ in Eq. 
(9), it was estimated in Ref. [22] by means of Padc approximations approach. 
The input information, used in these estimates, is the analytical result for the 
Uf ^-contribution to P4 calculated in [23] . In our normalization conditions it has 
the following form 



/5f 



1 



1024 



1205 
2916 



152 ■ 



-0.0017969nf 



The approximation of Ref. [22], which is taking it into account, reads: 
10' 

P4 



1024 
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(14) 

(15) 
5). 
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We will use this estimate in our studies. The related to the MS-scheme anoma- 
lous mass dimension function is defined as 



dlnmf, 



7m (Os) 
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(16) 
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The four-loop correction was calculated independently in Ref. [24] and in 
Ref. [25]. Let us present the explicitly known coefficients: 
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^ 98.933 - 19.108nf + 0.27616 nf + 0.005793 nf fa 11.034 (uf = 5) . 

In the same work [22] the following model for the five-loop coefficient of was 
proposed: 

74 fa 530 -143 nf + 6.67 nf + 0.037 nf- 8.54 X 10~^n| fa -13.68 (uf = 5) (21) 

It is based on application of the variant of the Pade approximation approach, 
used for getting Eq. (15) discussed above. The explicit analytical expression 
for the nf contribution to 74, extracted from the QED results of Ref. [26] and 
confirmed later on in [23], namely 



[4] 
74 



1/65 5C3 
12 V 5184 ^ 324 



^4 , - 
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(22) 



was used. Notice that the numerical value of Eq. (21) is negative. This means 
that the uncertainties of this Padc estimate are not small. Indeed, the analytical 
calculations of nf contribution to Eq. (21), performed in Ref. [27], gave 
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It does not agree with the similar coefficient in Eq. (21). Substituting Eq. (23) 
into Eq. (21) one can see that for Uf = 5 the estimate for 74 is still negative, 
but rather small, namely 74 ~ —4.76595. We will incorporate this value in our 
further analysis just by fixing some parts of existing five-loop ambiguities. We 
hope that this expression, obtained by matching the explicit results of Eq. (22) 
and Eq. (23) and the Pade resummation technique, may be improved in the 
future. 

3 Analytical continuation effects and APT approach 

Consider now in more detail the contributions to T^^^ from the Minkowskian 
coefficients, defined in Eqs. (3)-(5). Let us remind, that they are composed from 
the Euclidean terms (see Eqs. (6) (8)) and the combinations of the coefficients 
of the QCD RG functions [3{as) and 7m(as), which are defining proportional to 
TT^ analytical continuation effects in Eqs. (3)-(5). These kinematic effects turn 
out to be negative and quite sizable. Thus, the coefficients AFi {i > 2) are 
smaller, than their Euclidean "analogs" df . Indeed, for Uf = 5 we have [12] 

AFsRi 29.147. (24) 

The numerical values of other terms obey the same pattern [10]: 

ATa « 41.758; (25) 
Ar4 K. -825.75. (26) 

This means that it is of interest how the effects of analytical continuation influ- 
ence other results of perturbative QCD predictions. In the beginning of 80s this 
problem was discussed in the number of works on the subject (see, e.g., [28, 29]). 
At that time the problem of resummation of the 7r^-contributions to Tjjbb was 
also considered in Ref. [11]. However, the real interest to rcsummations of the 
analytical continuation effects was attracted later on after appearance of Contour 
Improved Technique (CIT) [30,31] and Analytic Perturbation Theory (APT), 
in particular. This method was proposed and developed by D. V. Shirkov and 
I. L. Solovtsov in the process of common investigations in Rcfs. [32]-[34] and in 
the separate publications as well (see works done with the decisive contribution 
of I. L. Solovtsov, [35, 36], and the ones created by inspiration of D.V. Shirkov 
[37,38]). This QCD approach was already used in various applications (see, 
e.g., [39,40]). Among them are the studies of Higgs boson decay into bb-pair 
with the help of Fractional Analytical Perturbation Theory (FAPT) [41], which 
are complementary to definite considerations of Ref. [42]. Quite recently some 
new [43] and even a bit corrected [44] discussions of applications of FAPT to 
Fjjbb s-iid their comparisons with the results of Ref. [42], appeared in the liter- 
ature. In view of this, we will focus ourselves here on the brief discussions of 
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several ideas of Rcf. [42]. It is worth to stress, that this work was motivated in 
part by the desire to understand whether the bridge may be built between the 
renormalon approach (for a review see, e.g., [45]) and the resummation of the 
proportional to /?o-efFects within Shirkov-Solovtsov APT. Note, that one of the 
main cornerstones of renormalon approach is /3o-resummation procedure as well. 
Today, thanks to the works of Refs. [41,44] it is understood, that this bridge 
does exist. The essential point in the studies of Ref. [42] is that the CIT [30] and 
the concept of b-quark invariant mass are playing an important role. In view of 
the fact that the positive features of the invariant mass arc sometimes not taken 
used to the total extent, let us remind the basic steps of definitions of this QCD 
parameter: 

1. Define the running quark mass through the solution of the following RG 
equation for the anomalous dimension term (AD): 



mb(MH) = mb(mb) exp 



™2/ 



- 2 



as{mb) 



dx 



(27) 



2. Take the integral in the r.h.s. of Eq. (27): 



m^(MH)=m^(mb) 
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4. Calculate its coefficients, expressed through 
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5. Define the invariant mass 



rhb = mb(mb) 



To. 



as(mb)''o^-D(as(mb)) 



(32) 



It should be stressed, that this definition seems to be simpler than one, 
introduced in Ref. [15], namely 



ihb = mb(mb) 



70 



(2/3oa,(mb))'5o^£)(a,(mb)) 



-1 



(33) 



which is often used in the literature. 



In the large-/?o approximation the expression for Tjjbb can be expressed as [42] 
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(34) 



(35) 



where 6n = n + uq — 1, vq = 270 //3o- Taking now into account that within 
large-/3o approximation, one has the following LO-expression 



as(MH)LO 



/?oln(M^/A2) 



(36) 



fixing n = and expanding Aq to first order in we arc getting the function: 



An 



1 sin(6 arctan(7r/LMH; 



(37) 



where h = - I, Lm^ = \n{M^/A^). This expression was first derived in 
80s in Ref. [11]. Unfortunately, its usefulness was not understood at this time. 
At the new stage of the development of QCD the analogs of this formula, are 
forming the basis of FAPT method [41] , which is allowing to resum not only the 
terms proportional to 70 and /Sq, but higher order corrections of RG functions 
as well. Thus, the extension of the Shirkov-Solovtsov method to the case of 
fractional powers [41] may be considered in part as the generalization for CIT 
resummation of "large-/?o" contributions, which also arises within renormalon 
approach [42]. In view of the appearance of the works of Refs. [41, 44], [43], it 
may be interesting to compare in the future the results of these two approaches 
in more detail. 
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4 On-shell and RG-resummation approach 



4.1 On-shell parameterization 

In the previous section we derived the relation between running and invariant b- 
quark masses. However, there is also the possibility to use on-shell approach and 
express the width T^^^ through the b-quark pole mass mt and the MS-scheme 
coupling constant Q!s(Mh) in different orders of perturbation theory and com- 
pare the results obtained with the running mass motivated RG-resummation 
approach. This analysis was done at the a^-level in Refs. [46,47]. In these 
studies the effects of mass dependent 0(asm^/M^)-corrections, extracted from 
the calculations of Ref. [48], were also included. Among most important results 
of Refs. [46, 47] were the explicit demonstration of the importance of taking into 
account a^-corrections in the on-shell approach. Indeed, these effects turned out 
to be rather important for decreasing the difference between the T^^^ expres- 
sions, evaluated in the on-shell and RG-rcsummed approaches. The results of 
our studies were confirmed later on by the considerations of Ref. [49] , where the 
expression of the 0(a^m^/M^)-contribution was evaluated and included. Note, 
however, that for the considered at present masses of Higgs boson these effects 
are less important, than higher order pcrturbativc QCD corrections, and can 
be safely neglected. Keeping in mind the demands of Tevatron and LHC ex- 
periments and the ongoing discussion of the scientific program for International 
Linear Collider, in this section we will study the similar problem in more detail, 
taking into account the information on available at present higher-order QCD 
corrections to the RG-functions, coefficient function for T^^th Eq. (1) (see 
Sec. 2) and the relation between running and on-shell masses which we present 
in the following form: 

mb(mb) _ -, _ 8 
ml "3 



- l--a5(mb) -18.556 a5(mb)^-175.76a5 (mb)^- 1892 as (mb)"^. (38) 
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The a^-correction is the result of calculations of Ref. [50], confirmed later on 
in [51]. The a^-term comes from the analytical calculations of Ref. [52] and is 
confirming semi-analytical similar result, obtained in Ref. [53]. Note, that its 
coefficient turned out to be in a good agreement with the ECH/PMS estimate 
of Ref. [14]. This fact and the success of the ECH/PMS prediction for the value 
of df term [10], wc arc supplementing the estimates of Ref. [14] by definite 
RG-inspired considerations and get our personal ECH-inspired number for the 
coefficient of a^-term in Eq. (38). Proceeding further on with the help of the 
derived in Ref. [14] RG equations for the transformation of m^(MH) to m^(mb) 
and of as (mb) to as(MH), we get the following analog of Eq. (1): 



r - — 
^ Hbb — '-o 



l+AT\as (Mh) + Ar^as (Mh)^ + Ar^a^ (Mh)^ + Ar^a^ (Mr) 



(39) 
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where Fq = 3V2/(87r) Gp Mr and 



3-2L; (40) 

Ar^ = -4.5202- 18.139 L + 0.08333 L^; (41) 

Ar| = -316.878 - 133.421 L - 1.15509 + 0.050926 ; (42) 
Ar^ = -4366. 17- 1094.62 L- 55.867 L^-l. 8065 + 0.04774 (43) 



and L = ln(Mg/m^). We will define now the QCD coupling constant in different 
orders of perturbation theory as 



as(MH)NLO 

a5(MH)N2LO 
a,(MH)N3LO 
a,(MH)N4LO 
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5 
3^ 



Pa 
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where Log = 2 ln(MH/A|^~^'') and the additional terms Aa5(MH)N3LO ^"^^ 
Aas(MH)N4LO were obtained in Refs. [54, 55] with the corresponding N^LO and 
N^LO matching conditions, which are allowing to determine the values of A^~^^ 

from A^~^^ by passing threshold of production of heavy flavor (in our case on- 
shell mass mb). The analytical results of Ref. [55] are in complete agreement with 

the mixture of previous analogous analytical and semi-analytical calculations of 
Refs. [56, 57]. These conditions generalize to higher orders the NLO and N^LO 
formulae, derived in Ref. [58] (the corresponding N^LO relation was corrected 
in Ref. [59]). To save the space, we will not present here the explicit form of 
these equations. An interested reader can consult Ref. [60], where the results 
of Refs. [58, 59] and [54] are presented for the case of considering b-quark pole 



10 



mass as the matching point. The corresponding N^LO version will be presented 
elsewhere. In order to perform the on-shell analysis we are using the results for 
the b-quark pole mass values from Ref. [61], which are increasing from LO to 
N^LO, and the LO, NLO, N^LO expressions for Ai^"^\ related to one of the 
most recent parton distribution fits Ref. [62]. It can be shown, that NLO and 
N^LO values of these parameters are in agreement with the central values of 
A^~^^ obtained in the process of the performed in Ref. [60] fits of Tevatron vN 

deep-inelastic scattering data. In view of this the N^LO value of A^"^'' will be 

^ MS 

taken from Ref. [60] . Leaving aside the discussions of the possible convergence 
of the fits, at the N^LO we will use the same value of A^"^-* as at the N^LO 
and will make the guess about the value of the b-quark mass. Using now the 
matching conditions of Ref. [55] transformed into the form, given in Ref. [60], 
we obtain the following values of A^"^'* : 



order 


mb GeV 


a£l-) MeV 


A^L-^ MeV 


LO 


4.74 


220 


168 


NLO 


4.86 


347 


254 


N^LO 


5.02 


331 


242 


N^LO 


5.23 


333 


243 


N^LO 


5.45 


333 


241 



Table 1: The values of the parameters on the on-shell studies. 

In order to study the effects of separate NLO, N^LO, N-'^LO and N^LO 
corrections to different parameterizations, it is also useful to define the quantity 

%bb = rHbb/ro (48) 

and study its Mh dependence. At the final step it is necessary to substitute the 
results from Table 1 into Eq. (39)-(47) by iterative way and compare necessary 
curves. They are shown in Fig. 1 both for Ljj^^, and -RnBb quantities. 

4.2 Renormalization group resummation approach 

In order to apply RG-resummation approach to Lg^b it is necessary to combine 
definition of Eq. (1), the basic formulae (27)-(29), the concrete numbers of 
Eqs. (30)-(31), the transformation relation of Eq. (38) and the expressions for 
as(MH) defined through Eqs. (44)-(47) and the similar expressions for as(mb). 
The results are presented at Fig. 2. 
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Fig. 1: The analyzed quantities in the OS-approach. 
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Fig. 2: The analyzed quantities in the RG-approach. 
Fjjgj, may be presented in the following form 



1 + j;AFi4(MH) 



i>l 



(49) 



X (1 - -as(mb) - 18.556 as (mb)^ - 175.76 a5(mb)^ - 1892.2 as(mb)'') 



where 



AD{asf = 1 + 2.351 + 4.383 + 3.873 - 15.153 a 



(50) 
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Different curves at Figures 1 and 2 arc related to applications of the results for 
Fjj^j, and -RfjEb °f ^^e step-by-step substitution into the definitions of the corre- 
sponding expressions for the QCD coupling constants a5(MH) and 05(0)3) from 
Eqs. (36), (44)-(47). The considerations of Fig. 1 confirm the findings of Refs. 
[46, 47] on the importance of the 1- and 2-loop contributions to Fjjgj^ in parame- 
terization with on-shell mass. On the other hand, the other two corrections are 
not so big, but they both have a tendency to decrease the value of T^^^ to its 
RG-improved expression (see Fig. 2). Moreover, the RG-resummation approach 
demonstrate clearly, that the perturbativc theory to this quantity are well un- 
der control. The r.h.s. parts of Fig. 1 and Fig. 2 demonstrate the stability of 
good convergence of the related perturbative approximations. Thus, the consid- 
eration of the order corrections, calculated in Ref. [10], support the feature 
of minimizing the difference between on-shell and RG-resummed parameteriza- 
tions, already observed at the a^-level in Rcf. [63]. The l.h.s. plots of Fig. 1 and 
Fig. 2 are more interesting from phenomenological point of view. It will be the 
next task to compare these results for the Higgs decay width, obtained by using 
FAPT approach of Ref. [41] and resummed FAPT analysis of Ref. [43]. Thus, our 
considerations leave space for further studies of peculiar features (presumably, 
better convergence) of modified APT predictions. 
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